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Abstract 

The structure of the second + state of 10 Be is investigated using a mi- 
croscopic a+a+n+n model based on the molecular-orbit (MO) model. The 
second + state, which has dominantly the (l/2 + ) 2 configuration, is shown 
to have a particularly enlarged a-a structure. The kinetic energy of the two 
valence neutrons occupying along the a-a axis is reduced remarkably due to 
the strong a clustering and, simultaneously, the spin-orbit interaction unex- 
pectedly plays important role to make the energy of this state much lower. 
The mixing of states with different spin structure is shown to be important in 
negative-parity states. The experimentally observed small-level spacing be- 
tween 1~ and 2~ (~300 keV) is found to be an evidence of this spin-mixing 
effect. 12 Be is also investigated using a+a+4n model, in which four valence 
neutrons are considered to occupy the (3/2") 2 (l/2 + ) 2 configuration. The 
energy surface of 12 Be is shown to exhibit similar characteristics, that the re- 
markable a clustering and the contribution of the spin-orbit interaction make 
the binding of the state with (3/2~) 2 (l/2 + ) 2 configuration properly stronger 
in comparison with the closed p-shell (3/2~) 2 (l/2 - ) 2 configuration. 

PACS number(s): 21.10.-k, 21.60.Gx 



I. INTRODUCTION 

Numerous experiments using unstable nuclear beams have succeeded to reveal exotic 
properties of /3-unstable nuclei, including a neutron halo structure Especially, the 

shift of the closed-shell structure is a characteristic behavior of systems with weakly bound 
neutrons. Recently, the contributions of the sd shell have been analyzed in N = 8 nuclei 
based on the shell model. A calculation has shown that the slow /3-decay of 12 Be to 12 B can 
be explained by an admixture of the sd shell in 12 Be (N = 8) in which the closed p-shell 
component must be less than 30% ||. This shows that the concept of magic numbers is 
vague in 12 Be. 

On the other hand, an a-a structure is well established in the Be and B region. Typi- 
cally in 9 Be, a microscopic a+a+n model has reproduced the properties of low-lying states 
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including the level inversion of l/2~ in the p-shell and l/2 + in the sd shell Also in 

10 Be, the microscopic a-cluster models have been applied 0,§, and we have quantitively 
shown the mechanism for the lowering of the sd orbits related with the clustering of the core 

According to our previous results for 10 Be, the main properties of this nucleus have been 
well described by the a+a+n+n model, where the orbits for the two valence neutrons are 
classified based on the molecular-orbit (MO) model. The dominant configuration of the two 
valence neutrons for the second + state is (l/2 + ) 2 , and the level inversion between l/2~ and 
l/2 + in 9 Be also holds in 10 Be. Here, the two valence neutrons stay along the a-a axis, and 
reduce the kinetic energy by enhancing the a-a distance (up to around 4 fm). Therefore, 
the l/2 + orbit in this case is not a spherical s-wave, but a polarized orbit with the c?-wave 
component. This feature of sd mixing in the l/2 + orbit can be qualitatively interpreted in 
terms of deformed models ([220] expression in the Nilsson diagram). 

The main purpose of this paper is to show the mechanism for the appearance of this 
(l/2 + ) 2 configuration in low-lying energy in detail. In our previous work, the calculated 
excitation energy of the second + state has been higher by ~5 MeV without the spin-orbit 
interaction than the experimental one 0. The spin-orbit interaction has been found to 
decrease the excitation energy by 3 MeV by adding a single optimal wave function with 
S — 1. Hence, it is necessary to analyze the contribution of the spin-orbit interaction more 
precisely to clarify the lowering mechanism of the second + state. The other reduction 
mechanism of the kinetic energy by more accurate description of the neutron-tail should 
also be taken into account. 

Firstly, the spin-orbit interaction is able to contribute to the second + state, when spin- 
triplet states for the valence neutrons are included among the basis states. If the l/2 + state 
consists of the pure s-orbit, naturally there is no contribution of the spin-orbit interaction. 
In Be isotopes, the state contains the (i-orbit component. However, the spin-orbit interaction 
again vanishes, when the two valence neutrons occupy along the a-a axis, since two neutrons 
with the same spatial distribution construct only the spin-singlet state. This has been the 
situation in traditional MO models f7]f|]. It will be shown that when one of the valence 
neutrons deviates from the a-a axis, the spin-triplet state can be constructed, and the 
spin-orbit interaction strongly acts between this state and the original (l/2 + ) 2 configuration 
with the spin-singlet. The extension of the model space also enables us to improve the 
description of the neutron-tail. We estimate the reduction of the kinetic energy due to the 
improvement of the neutron-tail. These improvements will be shown to largely decrease the 
excitation energy of the second + state, and that the calculated energy just corresponds to 
the experimental value. 

Using the present extended MO model, negative-parity states in 10 Be are also investi- 
gated. The mixing of different spin structure (S — 0, S — 1) is found to also be important 
in the negative-parity states. The i^-mixing effect strongly acts on the 2~ state, and the 
calculated small-level spacing between 1~ and 2~ (~300 keV) agrees with the experimental 
value. 

Also, in 12 Be these two important mechanism for lowering the (l/2 + ) 2 configuration 
are examined. In 12 Be, the configuration for four valence neutrons around the a-a core is 
considered to correspond to the two pair configurations, (3/2 - ) 2 and (l/2 + ) 2 . These two pair 
configurations appear in the ground + state and the second + state of 10 Be, respectively. 
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The energy surface for the ground state of 12 Be will be shown, in which the spin-orbit 
interaction is seen to make the binding for the (3/2~) 2 (l/2 + ) 2 configuration properly stronger 
in comparison with the (3/2~) 2 (l/2~) 2 configuration. This is related to the breaking of the 
N = 8 (closed p-shell) neutron-magic number. 

This paper is organized as follows. In Sec. II, we summarize a description of the single- 
particle orbits around the two a clusters based on the MO model. In Sec. Ill, the contri- 
bution of the spin-orbit interaction to the (l/2 + ) 2 configuration and the reduction of the 
kinetic energy due to an improvement of the neutron-tail are considered. These effects are 
discussed in detail for the second + state of 10 Be (A), for the negative parity states of 10 Be 
(B), and for 12 Be (C). The conclusion is given in Sec. IV. 



II. EXTENDED MOLECULAR ORBITAL MODEL 



We introduce a microscopic a+a+2n model for 10 Be and a+a+4n model for 12 Be. The 
neutron configurations are classified according to the molecular-orbit (MO) picture [pTOfl . The 
details of the framework are given in Ref. 0; here, the essential part is presented. 

The Hamiltonian is the same as in Ref. [[J, and consists of a kinetic-energy term, a 
central two-body interaction term, a two-body spin-orbit interaction term, and a Coulomb 
interaction term: 



h = E^-^ + E^' + E^ s + E^-( 1 -^)( 1 



i<j 



i<j 



i<j 



(1) 



The effective nucleon-nucleon interaction is Volkov No. 2 [11] for the central part and the 



G3RS spin-orbit term [|l!| for the spin-orbit part, as follows: 



Vij = {V x e- air ^ - V 2 e~ a2r ^}{W - MP a P T + BP a - HP T }, 



(2) 



V$ = Vtf*{e- air « - e- a < }L ■ SP 31 , (3) 

where P31 is a projection operator onto the triplet odd state. The parameters are V\ = 
-60.650 MeV, V 2 = 61.140 MeV, ^ = 0.980 fm~ 2 and a 2 = 0.309 fm" 2 for the central 
interaction, and Vq s = 2000 MeV, ai = 5.00 fm -2 , and a 2 = 2.778 fm -2 for the spin-orbit 
interaction. We employ the Majorana exchange parameter, M = 0.6 (W = 0.4), the Bartlett 
exchange parameter, B = 0.125, and the Heisenberg exchange parameter, H = 0.125, for 
the Volkov interaction (using B and H, there is no neutron-neutron bound state). All of 
these parameters are determined from the a + n and a + a scattering phase shifts and the 
binding energy of the deuteron || . 

The total wave function is fully antisymmetrized and expressed by a superposition of 
terms centered to different relative distances between the two a clusters (d) with various 
configurations of the valence neutrons. The projection to the eigen-states of angular mo- 
mentum (J) is performed numerically. All nucleons are described by Gaussians with the 
oscillator parameter (s) equal to 1.46 fm. The a clusters located at d/2 and — d/2 on the 
z-axis consist of four nucleons: 

(Q) = G&GtG&G&XptX*iXniXni. (4) 
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G represents Gaussians: 

GR a = (^) 4 exp[-i/(f - R a )\ v = l/2s\ R a = {de z /2, -de z /2}. (5) 

Each valence neutron ((p C iXa) around the a-a core is expressed by a linear combination of 
local Gaussians: 

<f>ciXci = J29jG R i i Xti, (6) 

3 

3 

G Ri = (^exph^f - Rl)% v = l/2s\ (7) 

The level structure is solved by superposing basis states with different a-a distances and 
configurations after the angular momentum projection. 

In the MO model, valence neutrons are expressed by a linear combination of orbits around 
two a clusters. The orbit around each a cluster is called the atomic orbit (AO). Here, we 
note that the valence neutrons and neutrons in the a clusters are identical particles, and 
antisymmetrization imposes the AO forbidden space. The lowest AO has one node and 
parity minus, which is the p-orbit. When a linear combination of two AOs (p-orbit) around 
the left-a cluster and the right-a cluster is summed up by the same sign, the resultant MO 
also has negative parity and one node. In this case, the MO is restricted to spread along an 
axis perpendicular to the a-a (z) axis, which is the so-called 7r-orbit. It cannot spread along 
the z-axis where the two a clusters are already located. If two p-orbits are summed up by 
different signs, the resultant MO has two nodes and positive parity. This MO can spread 
to all directions, and the optimal direction becomes the z-direction. This is the so-called 
cr-orbit, and the energy becomes lower as the distance between two a clusters increases. 

In Be isotopes, the low-lying orbits with K n = 3/2~, l/2 + , and l/2~ are important; 
these are classified based on the MO model. The l/2 + state is the cr-orbit, and the spin- 
orbit interaction splits the 7r-orbit to 3/2~ and 1/2 - . In the present framework, each valence 
neutron is introduced to have a definite K n at the zero limit of centers of local Gaussians 
{{R J n }) describing the spatial distribution of the orbit. The precise positions of {R 3 n } are 
determined variationally for each a-a distance before the angular- momentum projection. 
Since the values of {R^} are optimized to be finite, the orbits are not exactly the eigen-state 
of K 77 , and are labeled as K 7 ". 

In 10 Be, when two valence neutrons occupy these single particle orbits, three configura- 
tions with the total K = K\ + K 2 = are generated. For the ground state, valence neutrons 
occupy orbits with K n = 3/2~ and K 7 " = — 3/2 _ ; these 7r-orbits are expressed as linear 
combinations of p x and p y : 

$(3/2~, -3/2-) = A[<ptU ( 2 } (0 c iXci)(0 c2 X C 2)]. (8) 

0clXcl = {(Px + iPy)+a + (Px + ip y )-a}\n |), (9) 
<Pc2Xc2 = {(Px ~ iPy)+a + (Px ~ iPy)-a}\n |). (10) 
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Each MO is constructed from two AOs around a left and right a cluster. We introduce 
a variational parameter (a) on the ol-ol axis (z-axis); and (p x ± ip y ) +a is p x ± ip y , whose 
center is +a on the z-axis; and (p x ± ip y )- a is p x ± whose center is —a on the z-axis. 
Here, the spatial part and the spin part of K in Eqs. (9) and (10) are introduced to be 
parallel, for which the spin-orbit interaction acts attractively. The spin-up valence neutron 
has K\ = 3/2~ (rY"n|n f )), an d the spin-down valence neutron has K\ = 3/2~ (rYi_i|n |)). 
There (p x ) +a , (p x )- a , (p y ) +a , and {p y )- a are expressed as the combination of two Gaussians, 
whose centers are shifted by a variational parameter (b) perpendicular to the z-axis. For 
example, 

(Px)+a G aez+be x ^ae*2— be^; i.Py)+a G a e z +be y ^ae z — be y - (H) 

These parameters a and b, are optimized by using the Cooling Method in antisymmetrized 
molecular dynamics (AMD) [|T3"|-|i7| for each a-a distance. 

The a-orbit pair, (l/2 + ) 2 , is prepared, where each u-orbit is expressed by subtracting 
two AOs at +a = d/2 and —a = —d/2: 

$(l/2 + , -l/2 + ) = ^[0i a) 0? ) (0 c iXci)(toc2)] ) (12) 

0ciXci = {(p)+a ~ (p)-a}\n t) a = d/2, (13) 

C 2Xc2 = {(p)+ a - (p)- a }\n |> a = d/2, (14) 

where, 

(S±a = (^-G_5) ±B . (15) 

Since the a-orbit has two nodes, the direction of the orbit is not limited due to the Pauli 
principle. It can thus take all directions. Therefore, the direction and value of parameter b 
for the center of the Gaussians are optimized. As a result, it is optimized to the direction 
of the a-a axis. This basis state shares the dominant component of the second + state. 

III. RESULTS OBTAINED USING THE EXTENDED MO MODEL SPACE 

A. J T = 0+ states in 10 Be 

Before giving the results, we summarize the calculated energies and the optimal a-a 
distance for three configurations ((3/2 - ) 2 , (l/2~) 2 , and (l/2 + ) 2 ) discussed in Ref. [[J. In 
Fig. 1, the 0+ energy curves of $(3/2~, -3/2"), $(l/2 + , -l/2 + ), and $(l/2~, -1/2") in 
10 Be are presented as a function of the a-a distance (d). These configurations are the 
dominant components of the first, second, and third + states in 10 Be, respectively. 



Fig. 1 
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The energy of $(1/2 , — 1/2 ) (dominant component of 0^) becomes lower as d increases, 
and the energy becomes minimum at an a-a distance of 4 ~ 5 fm. However, it should be 
noted that the energy difference between $(l/2 + , — 1/2 + ) and $(3/2 _ , — 3/2~) in Fig. 1 is 
about 11 MeV. This is much larger than the experimental excitation energy of the second 
0+ state (6.26 MeV 0) by 5 MeV. 

There are two reasons for this lack of binding energy by 5 MeV in the second + state. 
One is that the contribution of the spin-orbit interaction is not taken into account; the 
other is that the tail of the valence neutron is not correctly described. As for the spin-orbit 
interaction, in the traditional MO model, the spin-orbit interaction has not contributed to 
the + state with this configuration. This is because, the two valence neutrons occupy the 
same spatial configurations with opposite spin direction, and construct a spin-singlet state. 
Since it is impossible for two valence neutrons with the same spin direction to occupy the 
same a-orbit along the a-a axis, the spin-triplet state is prepared by making one of the 
valence neutrons deviate from the a-a axis. As for the tail effect of the valence neutrons, in 
this analysis, the tail behavior is expressed by the superposition of local Gaussians. 

These two improvements are shown for the $(l/2 + , — 1/2 + ) configuration with the opti- 
mal a-a distance of 4 fm. The optimized parameter b is 2.02e 2 (fm) ||. As listed in Table 
I, the spin-singlet $(l/2 + , — 1/2 + ) configuration has a + energy of —46.3 MeV. As shown 
in Eqs. (12)-(15), each valence neutron is described as a linear combination of four local 
Gaussians. When we optimize their linear combination, the energy slightly decreases to 
-47.7 MeV (noted as S = (a) 2 in Table I). 



Table I 



Next, the contribution of the spin-orbit interaction is taken into account to this spin-singlet 
state by preparing the spin-triplet state (S z = 1). One of the valence neutrons (0 c iXci) 
occupies the cx-orbit, and another one (0 C 2Xc2) with the same spin-direction deviates from 
the a-a axis. The later one is described by a local Gaussian centered at R; 4> C 2Xc2 = Gg\n f). 

— * 

The calculated energy surface for the second + state as a function of the parameter R is 
given in Fig. 2. The orthogonality of the ground + state to the ground + state is ensured 
in this calculation. The point on the x-z plain in Fig. 2 shows the position of the parameter 
R, and the contour map shows the total energy obtained by taking into account this coupling 
with the spin-triplet states. 



Fig. 2 



The lowest energy of —50.9 MeV is seen at R = (x,z) = (1.0,2.0) (fm) on the energy 
surface. The contribution of the spin-orbit interaction decreases the energy by more than 
3 MeV in comparison with the case of the spin-singlet state. Although the energy of such 
a spin-triplet state, itself, is much higher than that of $(l/2 + , — 1/2 + ) by about 20 MeV, 
the large coupling between the spin-singlet and spin-triplet states enables the spin-orbit 
interaction to attractively act. As a result, more than half of the under-binding by 5 MeV 
for the second + state is explained by this spin-orbit contribution. 
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— * 

To see this effect more correctly, a lot of spin-triplet states with different R values are 
adopted. The employed R values are 0.0, 1.0, 2.0 (fm) for the x-direction (perpendicular 
to the a-a direction) and 0.0, 1.0, 2.0, 3.0, 4.0 (fm) for the z-direction (the a-a direction). 
In Table I, the column noted as spin-orbit gives the energy of the second + state when 

— * 

these 15 spin-triplet states with different R values are superposed. The calculated energy is 
—51.2 MeV; this is lower than the energy minimum point in Fig. 2 by only 300 keV. This is 
because, the spin-triplet basis states around the core have a large overlap with each other, 
and those apart from the core do not contribute to this coupling. This result shows that the 
contribution of the spin-orbit interaction for the (l/2 + ) 2 state is approximately taken into 
account when we employ only one Slater determinant on the energy surface. 

In this calculation, many spin-triplet states are employed; however, it should be noted 
that the spin-singlet state is only represented by the $(l/2 + , — l/2 + ) configuration. Thus, 
to improve the tail of the valence neutron in the second + state of 10 Be, it is necessary 
to superpose further the S z = basis states with many different R values, which contain 
the spin-singlet component. The distance of a-a is 4 fm, where one of the valence neutrons 
(0ciXei) occupies the a-orbit, and the other valence neutrons (0 C 2Xc2) with the opposite 
spin-direction is described by local Gaussian centered at R. The later one has the form of 
4>c2Xc2 — G^\n I). We superpose 15 basis states with S z = 0, and the employed R values 
are 0.0, 1.0, 2.0 (fm) for the x-direction (perpendicular to the a-a direction) and 0.0, 1.0, 
2.0, 3.0, 4.0 (fm) for the z-direction (the a-a direction). As listed in Table I, the energy 
becomes —52.7 MeV. These basis states increase the binding of the second + state by 1.5 
MeV compared with —51.2 MeV for the third row. Here, it is noted that these 5*2 = 
basis states contain both the spin-single and spin-triplet components. However, the increase 
of the binding energy by adding these S z — basis states is due to the spin-singlet basis 
states, since it is considered that the functional space of the spin-triplet states are already 
exhausted by the S z = 1 states. These results show that by taking into account both the 
spin-orbit interaction (spin-triplet basis states) and the tail effect (spin-singlet basis states), 
the under-binding for the second + state by 5 MeV is fully explained. The calculated 
excitation energy of the second + state agrees with the experimental one. 

To confirm the component of the spin-triplet configurations introduced, the expectation 
value of the spin-square (S 2 ) in the second + state of 10 Be was calculated. The value is 
0.26. This means that the spin-triplet configuration mixes in this state by about 13%. 

In this analysis, the a-a distance is restricted to 4 fm, but we now examine the contri- 
bution of the spin-orbit interaction with respect to the a-a distance. 



Table II 



Table II gives the contribution of spin-orbit interaction for the second + state with the a-a 
distance of 3 fm, 4 fm, and 5 fm. The column noted as S = (a) 2 gives the energy for 
the $(l/2 + , — 1/2 + ) configuration, where the linear combination of Gaussians is optimized. 
The column noted as +S Z = 1 gives the energy when we include the spin-orbit interaction 
by employing the spin-triplet basis states. The results show that the smaller is the a-a 
distance, the larger is the contribution of the coupling with the spin-triplet state. When the 
a-a distance is 5 fm, the coupling increases the binding energy by about 3 MeV; however, 
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when the a-a distance is 3 fm, the coupling increases the binding energy by about 4.5 MeV. 
Therefore, the coupling with the spin-triplet becomes stronger as the a-a distance becomes 
small. The column noted as +S Z = 0, 1 gives the energy when we further include the spin- 
singlet basis state to express the tail of the valence neutron. This effect does not have a 
strong dependence on the a-a distance, and almost constantly decreases the binding energy 
by 1.5 MeV. 

Finally, the calculation where all of these basis states with different a-a distance are 
superposed is performed based on the generator coordinate method (GCM). When we su- 
perpose all of the basis states in Table II, the calculated energy of the second + state 
becomes —54.3 MeV. Here, we superpose the basis states with a-a distances of 3, 4, 5, and 
6 fm, and the spin-triplet states are represented by one Slater determinant for each a-a 
distance, which correspond to the minimum point on the energy surface (Fig. 2). This cal- 
culated energy of —54.3 MeV is lower than a result calculated with only the <3>(l/2 + , — 1/2 + ) 
configuration (employed a-a distances are the same) by 4.6 MeV. 

B. Negative-parity states in 10 Be 

We examine this spin-mixing effect for the negative-parity states in 10 Be. Here, one of 
the valence neutrons is introduced to occupy a 7r-orbit with negative parity; another one 
occupies a <r-orbit with positive parity. In this way, two configurations with S z — and 
S z = 1 can be constructed. The S z — basis state ($(3/2~, — 1/2 + )) has K w = \~ . 

$(3/2~, -l/2 + ) = ^[0S Q) 0^(0 c iXci)(0 C 2X c2 )]. (16) 
c iXci = {(Px + ip y )+a + (Px + ip y )- a }\n |>, (17) 

<Pc2Xc2 = {(p) + a-(p)-a}\ni). (18) 

This basis state gives —47.7 MeV for 1~ and —46.7 MeV for 2~ at the optimal a-a distance 
of 3 fm. This energy splitting between 1~ and 2~ by 1 MeV is much larger than the 
experimental value of 300 keV (experimentally, the excitation energies of the \~ state and 
the 2~ state are 5.96 MeV and 6.26 MeV, respectively). At first, we only examine the tail 

— * 

effect of the valence neutrons with S z = 0. We introduce a parameter R for <f) C 2Xc2, which 
originally occupies the a-orbit: 

4>c2Xc2 = G n \n[). (19) 

However, even if we superpose states with many different positions (R), the energy splitting 
between these states is still about 1 MeV (1~ state is —49.9 MeV and the 2~ state is —48.9 
MeV). The employed R values are 0.0, 1.0, 2.0 (fm) for the x-direction (perpendicular to 
the a-a direction) and 0.0, 1.0, 2.0, 3.0, 4.0 (fm) for the z-direction (a-a direction). 

Therefore, to improve this discrepancy, we next include the basis states with S z = 1; 
$(3/2~, l/2 + ). The S z — 1 basis state generates K n = 2~ band. As a result of the 
calculation, the energy of the 2~ state with K n = 2~ (spin triplet) is found to be accidentally 



8 



close to the energy of the 2~ state with K n = 1~ projected from the S z = basis states. 
Thus, the coupling effect between the S z = basis states and the S z = 1 basis states is larger 
for the 2~ state than for the 1~ state, and the energy splitting becomes comparable to the 
experimental one. To see this effect, we superpose many S z = 1 basis states characterized 

— * " — * 

by the Gaussian center R; 4> C 2Xc2 = G^\n j). When we include many states with different R 
values in our bases states, the energies become —51.6 MeV (1~) and —51.2 MeV (2 _ ). The 
employed R values are 0.1, 1.0, 2.0 (fm) for the x-direction and 0.0, 1.0, 2.0, 3.0, 4.0 (fm) 
for the ^-direction. The experimental small level splitting between 1~ and 2~ is found to be 
evidence of the spin vibration in the 2~ state. 



C. The structure of Be 



The large contribution of this spin-orbit interaction will be discussed concerning 12 Be. 
In the previous subsection, we have discussed for 10 Be that this effect is more important as 
the ol-ol becomes smaller, and in 12 Be, we show that the distance is smaller than that of the 
second + state of 10 Be. In 12 Be, four valence neutrons rotate around two a clusters and, 
mainly, two configurations are important for the + ground state One is (3/2 - ) 2 (l/2~) 2 
for the four valence neutrons, which corresponds to the closed p-shell configuration of the 
neutrons at the zero limit of the ot-a distance. The other configuration is (3/2~) 2 (l/2 + ) 2 , 
where two of the four valence neutrons occupy the cr-orbit. It is necessary to compare the 
energy of these two configurations as a function of the ol-ol distance. The configurations of 
(3/2-) 2 (l/2-) 2 and (3/2-) 2 (l/2+) 2 are constructed as follows: 

$(3/2-, -3/2-, -1/2-, 1/2") = ^[^ a) ^ a) (0ciX c i)(0 C 2Xc2)(0 C 3Xc3)(0c4X C 4)] (20) 



<\>c\Xc\ = {(Px + iPy)+a + (Px + iPy)- a }\n |), (21) 

<t>c2Xc2 = {(Px ~ iPy)+a + (Px ~ Wy)-a}\n |), (22) 

0c3Xc3 = {(Px ~ iPy)+a + (Px - ip y )-a}\n f), (23) 

</Wc4 = {(Px + iPy)+a + (Px + iPy)-a}\n I), (24) 



$(3/2-, -3/2-, 1/2+, -1/2+) = ^[0^ ) 0? ) (0clX C l)(0 C 2Xc 2 )(0c3Xc3)(0c4X C 4)] (25) 

4>clXcl = {{Px + iPy)+a + {Px + ip y )-a}\n T}' (26) 

0c2Xc2 = {(Px - iPy)+a + (Px ~ ip y )-a}\n |), (27) 

toes = {(p) +a - (p)- a }\n T) a = d/2, (28) 

<Pc4Xc4 = {(p) +a -(p)-a}\nl) a = d/2. (29) 

The dotted line in Fig. 3 shows + energies of (3/2-) 2 (l/2 + ) 2 and the dashed line shows 
(3/2-) 2 (l/2-) 2 with respect to the a-a distance. 
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Fig. 3 



When the a-a distance is small, for example 2 fm, the dominant configuration of the four va- 
lence neutrons is (3/2~) 2 (l/2~) 2 for the ground state, which corresponds to the closed p-shell 
configuration at the limit of the a-a distance, zero. On the other hand, the (3/2 _ ) 2 (l/2 + ) 2 
configuration for the four valence neutrons becomes lower as the a-a distance is increased. 
However, it is still higher than the (3/2~) 2 (l/2~) 2 configuration by a few MeV. Here, we 
show the importance of the coupling between the (3/2~) 2 (l/2 + ) 2 configuration and the 
spin-triplet state for the last two neutrons ((l/2 + ) 2 ). The spin-triplet state for (l/2 + ) 2 is 
introduced as follows: 

$(3/2-, -3/2-, l/2 + ,G M \ T» =^[0S a) 0? ) (0 c iXci)(0c2X c2 )(0 C 3Xc3)(0 C 4X C 4)] (30) 
0ciXci = {(Px + ip y )+a + (Px + ip y )- a }\n T), (31) 

0c2Xc2 = {(Px~ iPy)+a + (Px ~ Wy)-a}\n I), (32) 

<tac3 = {(p)+a - (p)-a}\n T) a = d/2, (33) 

C 4Xc4 = Gg\n t). (34) 
The coupling between the (3/2-) 2 (l/2 + ) 2 configuration and the spin-triplet states is shown 

— * 

in Fig. 4 as a function of the parameter R. 



Fig. 4 



The point on the x-z plain in Fig. 4 shows the Gaussian-center R in the spin-triplet state. 
The contour map shows the energy calculated by taking into account this coupling with the 
spin-triplet state. This a-a distance is 3 fm, and this optimal distance is smaller than that 
for the second + state of 10 Be by about 1 fm, since two of the four valence neutrons occupy 
the (3/2-) 2 configuration and increase the binding energy of the system. Because of this 
effect, the coupling effect of the spin-triplet states is more important than in the case of the 
second + state of 10 Be. The minimal point of the surface shows that the energy gain due 
to the couplings is about 4 MeV in 12 Be. 

The solid line in Fig. 3 shows the + energy calculated by taking into account the 
coupling effect between the (3/2-) 2 (l/2 + ) 2 configuration and the S z — 0, 1 basis states for 
the last two valence neutrons. Due to the spin-orbit coupling, the energy is almost the same 
as that of (3/2~) 2 (l/2~) 2 corresponding to the closed p-shell configuration. Furthermore, 
the energy of (3/2-) 2 (l/2 + ) 2 is suggested to become even lower than (3/2~) 2 (l/2-) 2 when 
the pairing effect between (3/2 - ) 2 and (l/2~) 2 is taken into account. These effects plays 
crucial roles in accounting for breaking of the N = 8 magic number. 
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IV. CONCLUSION 



The second + state of 10 Be has been shown to be characterized by the cx-orbit of the two 
valence neutrons in terms of the molecular orbit (MO) model. The two valence neutrons 
stay along the a-a axis (the l/2 + orbit) and reduce the kinetic energy by enlarging the a-a 
distance. This simple description for the second + state gives a higher excitation energy 
by 5 MeV compared to the experimental one. To improve the description, spin-triplet basis 
states which have not been included in the traditional MO models are prepared by allowing 
a deviation of the neutron orbit from one just along the a-a axis. Thus, the spin-orbit 
interaction can be taken into account, and the calculated second + state becomes lower by 
3.5 MeV. The precise description of the neutron-tail also decreases the excitation energy by 
1.5 MeV. The discrepancy between the experimental excitation energy and the calculated 
one is compensated by these two effects. 

The spin-mixing effects were studied for the negative parity state of 10 Be. If we restrict 
ourselves to the S z — basis state, this energy splitting between the 1~ state and 2~ state 
is about 1 MeV, which is much larger than the experimental value of 303 keV. As a result 
of adding basis states with S z — 1, where the spin directions of the two valence neutrons are 
the same, the Admixing occurs especially for the 2 _ state. The contribution of the S z = 1 
state is larger for the 2~ state. The resultant energy splitting becomes comparable to the 
experimental one. Therefore, the experimental small level splitting between 1~ and 2" is 
considered to be a result of the spin vibration induced by the spin-orbit interaction in the 
2~ state. 

The coupling with the spin-triplet basis states is also important in the case of 12 Be. 
Without the spin-triplet basis state, the energy of the configuration (3/2~) 2 (l/2 + ) 2 is much 
higher than that of the closed p-shall configuration ((3/2~) 2 (l/2~) 2 ) by 4 MeV. However, 
the energy of (3/2~) 2 (l/2 + ) 2 is drastically decreased by coupling with the spin-triplet states. 
This is because the effect becomes stronger as the a-a distance becomes shorter, and 12 Be 
has an optimal a-a distance around 3 fm, which is smaller than the second + state of 
10 Be by 1 fm. The study shows that an energy of (3/2 _ ) 2 (l/2 + ) 2 is almost the same as 
(3/2~) 2 (l/2~) 2 , or even lower. This effect is suggested to play a crucial role in accounting 
for the dissipation of the N = 8 magic number in 12 Be. It is an interesting subject to 
analyze the binding mechanism and properties of the ground state by taking into account 
the pairing mixing among states with configurations of (3/2~) 2 (l/2~) 2 , (3/2~) 2 (l/2 + ) 2 , and 
(l/2~) 2 (l/2 + ) 2 . A detail analysis is going to be performed not only for this state, but also 
for excited states where new states with cluster structure have been recently observed. 
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FIGURES 



FIG. 1. The 0+ energy curve for $(3/2~, -3/2"), $(l/2~, -1/2"), and $(l/2 + , -1/2+) as a 
function of the a-a distance (d). 

FIG. 2. The energy surface of the second + state of 10 Be. The total energy including the 
coupling effect between the (l/2 + ) 2 configuration and the spin-triplet configuration is shown. The 
a-a distance is optimal 4 fm, and the spheres show a clusters. The axes (fm) show the position 
of the Gaussian-center of one valence neutron in the spin-triplet state which deviates from the 
cr-orbit. 

FIG. 3. The energy surface of the ground state of 12 Be. The total energy including the 
coupling effect between the (3/2~) 2 (l/2 + ) 2 configuration and the spin-triplet configuration for the 
last two valence neutrons is shown. The a-a distance is optimal 3 fm, and the spheres show a 
clusters. The axes (fm) show the position of the Gaussian-center of one valence neutron in the 
spin-triplet state which deviates from the cr-orbit 

FIG. 4. The 0+ energies for (3/2~) 2 (l/2 + ) 2 configuration (dotted line) and (3/2~) 2 (l/2-) 2 
configuration (dashed line) in 12 Be. In the solid line, the couplings among the (3/2~) 2 (l/2 + ) 2 
configuration, an optimal S z = 1 state, and S z = states are included. The resultant state almost 
degenerates with the (3/2~) 2 (l/2~) 2 configuration. 
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TABLES 



TABLE I. The energy convergence of the second 0+ of 10 Be. The a-a distance is fixed to 4 
fm. The row noted as 3>(l/2 + , — 1/2 + ) shows the energy of the 3>(l/2 + , — 1/2 + ) configuration, and 
the row S = (a) 2 shows the energy when a linear combination of local Gaussians is optimized. 
The row spin-orbit shows the energy including the coupling with the spin-triplet state. The row 
tail shows the calculated energy when the tail of the valence neutron is further expressed by the 
S z = basis states. The values in parentheses give the energy differecne with the previous row. 





number of Gaussians for one valence neutron 


+ energy 


$(l/2 + ,-l/2 + ) 


4 


-46.3 


S = (a) 2 


4 


-47.7 (-1.4) 


spin-orbit 


24 


-51.2 (-3.5) 


tail 


44 


-52.9 (-1.7) 



TABLE II. The second + energy of 10 Be for each a-a distance. The row S = (a) 2 shows 
the energy of the <5(l/2 + , — 1/2 + ) configuration, where a linear combination of local Gaussians is 
optimized. The row +S Z = 1 represents the energy including the coupling between 5 = (a) 2 and 
the spin-triplet state. The row +S Z = 0, 1 represents the energy including the coupling between 
S = (a) 2 and both the spin-singlet and the spin-triplet states. The values in parentheses give 
the energy differecne with the S = state. 



a-a distance (fm) 


S = {a) 2 (MeV) 


+S Z = 1 (MeV) 


+S Z = 0, 1 (MeV) 


3 


-44.5 


-48.9 (-4.4) 


-50.6 (-6.1) 


4 


-47.7 


-51.2 (-3.5) 


-52.9 (-5.0) 


5 


-46.9 


-50.0 (-3.1) 


-51.5 (-4.6) 
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6 1 — i — i — i — | — i — i — i — i — | — i — i — i — i — 

2 3 4 5 



alpha-alpha distance (fm) 



-50 



-51 - 




-6 n 1 1 1 1 1 1 1 1 1 

2 3 4 

alpha-alpha distance (fm) 



